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Abstract 

We compute the tree level cross section for gluon-gluon elastic scattering taking into account a 
dynamical gluon mass, and show that this mass scale is a natural regulator for this subprocess cross 
section. Using an eikonal approach in order to examine the relationship between this gluon-gluon 
scattering and the elastic pp and pp channels, we found that the dynamical gluon mass is of the 
same order of magnitude as the ad hoc infrared mass scale mg underlying eikonalized QCD-inspired 
models. We argue that this correspondence is not an accidental result, and that this dynamical 
scale indeed represents the onset of nonperturbative contributions to the elastic hadron-hadron 
scattering. We apply the eikonal model with a dynamical infrared mass scale to obtain predictions 
for cTtot^-, pPP'PP, slope BPP'PP, and differential elastic scattering cross section da^^/dt at Tevatron 
and CERN-LHC energies. 
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I. INTRODUCTION 



The increase of hadron-hadron total cross sections was theoretically predicted many years 
ago ^ and this prediction has been accurately verified by experiment Q]- At present the 
main theoretical approaches to explain this behavior are the Regge pole model and the 
QCD-inspired models. 

In the Regge pole model the increase of the total cross section is attributed to the exchange 
of a colorless state having the quantum numbers of the vacuum: the Pomeron j3j. In the 
QCD framework the Pomeron can be understood as the exchange of at least two gluons in a 

color singlet state Ol. A simple and interesting model for the Pomeron has been put forward 

n 

where it is evidenced the importance of the QCD nonperturbative vacuum 5]. One of the 
aspects of this nonperturbative physics appears in an infrared (IR) gluon mass scale which 
regulates the divergent behavior of the Pomeron exchange. 

In the QCD-inspired (or "mini-jet") models the increase of the total cross sections is 
associated with semihard scatterings of partons in the hadrons. The energy dependence of 
the cross sections is driven especially by gluon-gluon scattering processes, where the behavior 
of the gluon distribution function at small x exhibits the power law g{x,Q'^) ~ (see 



and references therein). In this case it is the gluon-gluon subprocess cross section 
that is potentially divergent at small transferred momenta. The procedure to regulate this 
behavior is the introduction of a purely ad hoc mass scale which separates the perturbative 
from the nonperturbative QCD region P, Q| . This mass scale, as well as the fixed coupling 
constant present in the elementary cross sections, are adjusted in order to obtain the best 
fits to the experimental data. 

On the other hand, several recent works have shown that the gluon may develop a dynam- 



ical mass (see the review jll| and the earlier work of Ref. [12|). This dynamical gluon mass 
was already successfully introduced in the Pomeron model of Landshoff and Nachtmann . 
Hence it is natural to ask if the arbitrary mass scale that appears in the QCD-inspired mod- 
els could be explained at a deeper level in terms of the dynamical gluon mass. This relation 
seems to be a plausible possibility and in this paper we will show that the dynamical gluon 
mass, as well as the IR finite coupling constant associated to it ^|, are in fact the natural 
regulators for the cross sections calculations. Since the behavior of the running coupling 
constant is constrained by the value of the dynamical gluon mass HQ , we will be able to 
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substitute the two ad hoc parameters in the "mini-jet" models p, |7|,|8|, namely the infrared 
mass scale (mo) and the effective value of the running coupling constant (os), by a physically 
well motivated one. In this way, beyond the natural interpretation of the arbitrary infrared 
mass scale in terms of a dynamical one, it is possible to decrease the number of parameters 
in this line of models. 

The paper is organized as follows: In the next section we introduce a dynamical gluon 
mass in the gluon-gluon scattering and compare the result to the standard one used in QCD- 
inspired models. In Sec. Ill we develop a QCD-inspired eikonal model in the light of the 
calculation of Sec. II. Our results are presented in the Sec. IV, where the best value for the 
dynamical gluon mass is determined and used thereafter to determine several quantities of 
pp and pp scattering. In Sec. V we present our conclusions. 



II. INFRARED MASS SCALE AND GLUON-GLUON ELASTIC SCATTERING 

In recent years there have been discussions in the literature about how to merge in a 
doubtless way the nonperturbative QCD results with the perturbative expansion. It is 
worth mentioning that Brodsky has several times called attention about the possibility to 
build up a skeleton expansion where the nonperturbative information would be included in 
vertices and propagators. In particular, the freezing of the QCD running coupling constant 
at low energy scales could allow to capture at an inclusive level the nonperturbative effects 
in a reliable way (see, for instance, Ref. (l5|). The freezing of the coupling constant and 
the existence of a dynamical gluon mass are intimately connected therefore they should 
appear systematically in this sought expansion. 

It is possible that such skeleton expansion could appear with the use of the pinch technique 
jioj l. With this technique the nonperturbative behavior of "gauge invariant" propagators 
and vertices could be computed nonperturbatively at one given order and substituted into 
the perturbative skeleton expansion. The fact that the "pinch" parts help to form gauge 
invariant quantities would result in well behaved matrix elements for the desired expansion. 

It is clear that we are still far from the kind of expansion discussed above, and we 
have to rely on more phenomenological approaches to go forward in this direction. One 
attempt to understand the effect of dynamically massive gluons was performed by Forshaw, 
Papavassiliou and Parrinello where they do introduce bare massive gluons and study 
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the amplitude behavior for some tree and one-loop level diagrams that could be relevant 
for diffractive scattering. In this approach, for example, the amplitude for the tree level 
process qq —* gg comes out with a mass dependence which is washed out in the high energy 
limit, but the massless limit is not recovered due to the presence of a numerically small 
mass independent term. The calculation is instructive but does not reproduce the high 
energy limit of massless gluons with two degrees of freedom, as discussed by Slavnov many 
years ago Actually the dynamical masses go to zero at large momenta and we should 
expect to recover the elementary cross sections of perturbative QCD in the high energy 
limit. Following this thought we could say that the sum over the polarizations should be 
performed as if the gluons were massless otherwise we would not map (at high energies) the 
desired skeleton expansion into the perturbative QCD expansion. 

According to the above discussion, we cannot work with a massive Yang-Mills theory or 
use a massive model where the third polarization state is provided by a massless scalar field. 
Thus, we will just assume a phenomenological procedure stated many years ago by Pagels 
and Stokar and named dynamical perturbation theory (DPT) 19]. The DPT approximation 
can be described as follows: amplitudes that do not vanish to all orders of perturbation 
theory are given by their free-field values. On the other hand, amplitudes that vanish in all 
orders in perturbation theory as oc exp {—1/g'^) {g is the coupling constant) are retained at 
lowest order. In our case this means that the effects of the dynamical gluon mass in the 
propagators and vertices will be retained, and the sum of polarizations will be performed 
for massless (free-field) gluons, because its signal (for massive gluons) will not vanish from 
the elementary cross section. In this approach, the differential elastic cross section for the 
process gg gg is written as 



da"'''",. 9to= 



s[AM^ -s-i] si i[AM^ - s - £] 

3 - 
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where and Mg are the expressions for the nonperturbative running coupling constant and 



for the dynamical gluon mass, respectively. They were obtained by Cornwall by means 
of the pinch technique in order to derive a gauge invariant Schwinger- Dyson equation for 
the gluon propagator. These expressions are given by 

"^^^'^ = /3oln [(g2 + 4M2(g2))/A2]' (2) 
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where Po = 11 — ^Uf {rif is the number of flavors) and A(= lyQco) is the QCD scale 
parameter. The latter expression has been determined as a fit to the numerical solution for 



the gluonic Schwinger- Dyson equation in the case of pure gauge QCD 



121. We also assume 



that the introduction of fermions does not change drastically this behavior. The gluon mass 
scale TUg has to be found phenomenologically, and a typical value is mg = 500 ± 200 MeV 
(for A = 300 MeV) [12, |20|. Note that we present the full gg gg cross section just for 
completeness, where the terms suppressed by powers of s are not important compared to 
leading In s perturbative corrections. However it should be pointed out that up to now 
higher order corrections have not been introduced in these type of models, and this is even 
more complicated if we consider the nonperturbative effects that we are introducing in this 
work. 

A different expression for the dynamical gluon mass can be found in Ref. [2^, given by 

M^il") = (4) 
which is consistent with the asymptotic behavior of Mg{q^) in the presence of the gluon 



condensates j2^. However, the calculation of the hadronic cross section does not depend 
strongly on the specific form of Mg(g^), but more on its IR value (i.e., the value of mg). 

In the limit ^ h? , the dynamical mass Mg{q'^) vanishes, and the nonperturbative QCD 
running coupling ots matches with the one-loop perturbative QCD one. Thus, in the limit 
of large enough g^, the expression (H)) reproduces its perturbative QCD counterpart: 



sis + t) St tis + t 



(5) 



To compute (0) we have used a vertex having a momentum dependent running coupling 
and a massive gluon propagator in the Feynman gauge, where the sum over gluon polariza- 
tions was performed for massless gluons. 

The total cross section a(s) = f-"^'^'' (da/di) di for the subprocess gg gg, that will be 
used in the next section to compose the eikonal term Xggi is obtained by integrating over 
— s < t < 0. In setting these kinematical limits we have neglected the momentum 



behavior in Eq. Q, as expected from our discussion on the weak dependence of hadronic 



cross sections on the specific form of M^(g ). A straightforward calculation yields 



2 



12s^ — 55mls^ + 12mis^ + 66m^s — 8m^ f s — 3m„ 

9 g 9 ^ _ 3 In ' 3 



Aim?gS[s — m^]^ \ m?g 



.2\ H 



(6) 



The asymptotic energy (s) dependence of the total cross section a^^'^{s) is of the following 
form 

We notice that the above result is similar to the asymptotic expression for the gluon-gluon 
total elastic cross section usually adopted in QCD- inspired models (QIM), 

^ S,, = (8) 

where the parameters mo and are assumed to be equal to 0.6 GeV and 0.5, respectively 
IliSl- We particularly call attention to these values, because they are of the same order of 
magnitude as the dynamical gluon mass scale (m^) and its frozen IR value of the coupling 
constant, obtained in other calculations of strongly interacting processes j2^. Therefore, 
all the point in here is how to connect these nonperturbative results to the straightforward 
perturbative QCD calculations. 



III. DYNAMICAL GLUON MASS AND QCD-INSPIRED EIKONAL MODELS 

A consistent calculation of high-energy hadron-hadron cross sections must be compatible 
with analyticity and unitarity constraints. The latter can be automatically satisfied by use 
of an eikonalized treatment of the semihard parton processes. In an eikonal representation, 
the total, elastic and inelastic cross sections are given by 

atot{s)=A'K I M6[l-e-^^(''^)cosxJ&,s)], (9) 
aei{s) = 2tt bdb\l- e-Xi(M+»x«(M)|2^ (10) 

aUs) = atotis) - cr,,{s) = 2n [ bdb[l- 6-'^'^''% (11) 

Jo 
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respectively, where s is the square of the total center-of-mass energy and x{b, s) is a complex 
eikonal function: s) = Xni^, s) + iXi{b, s). In this formalism, the factor e~'^^i^^'^^ in the 
expression (lllj) is interpreted as the probability that neither nucleon is broken up in a 
collision at impact parameter b. The ratio p of the real to the imaginary part of the forward 
scattering amplitude is given by 

Re{tjbdb[l-e'^('''^^]} 

Im{«/6c/6[l-e*x{M)]}' ^^^^ 
whereas the nuclear slope B and the differential elastic scattering cross section are given by 

lldbjl-e^ 



and 



= — [ bdb[l-e'^^'''^]Mqb)'' 



(14) 



dt ' ' 2n 

respectively, where Jo{x) is the Bessel function of the first kind. The eikonal function can 
be written as a combination of an even and odd eikonal terms related by crossing symmetry. 
In terms of the proton-proton {pp) and antiproton-proton (pp) scatterings, this combination 
reads x%ib, s) = s) ± s). 

Following the work of Block et al. [8|, we write the even eikonal as the sum of gluon-gluon, 
quark-gluon, and quark-quark contributions: 

X^{b,s) = Xqq{b, S) + Xqg{b, S) + Xgg{b, S) 

= l[aqg{s)W{b; /i,q) + a,g{s)W{b; Pgg) + agg{s)W{b; figg)]. (15) 

Here W{b]fi) is the overlap function at impact parameter space and o"jj(s) is the ele- 
mentary subprocess cross section of colliding quarks and gluons = q,g). The overlap 
function is usually associated with the Fourier transform of a dipole form factor, 

W{b;p) = ^{pbfKs{pb), (16) 

where K^lx) is the modified Bessel function of second kind. The W{b;fi) function is nor- 
malized so that J d'^bW{b; fi) = 1. The odd eikonal x~{b, s), that accounts for the difference 
between pp and pp channels, is parametrized as 

X-{b, s) = Com T.gg ^ e^'^ W{b- fi,,,), (17) 



where Egg is given by the expression (jHI) and mo is an arbitrary IR mass scale. Codd and 
fiodd are fitting parameters. We borrow this term, with its correct analyticity property, and 
write our odd eikonal as 

X-(6,s) = C-S^e-/^iy(6;/x-), (18) 

where rUg is the dynamical gluon mass and the parameters C~ and fi~ are constants to be 
fitted. The factor S is defined as 

E = H^£M, (19) 

which is just the expression (|7j) deprived of any momentum dependence, with the coupling 
constant ag set at its frozen IR value. This definition of S, when compared with the S^^ 
one, reveals explicitly the natural relation between the infrared mass scales mo and rrig. 

In the original Block et al. model the eikonal functions Xqqip^ s) and Xqgip, s), needed to 
describe the lower-energy forward data, are parametrized with terms dictated by the Regge 
phenomenology. Similarly, we parametrize our quark-quark and quark-gluon contributions 
as 

Tfi 

X,,{h,s)=i^A^W{h-^igg), (20) 

V ^ 



A' + B' In 



s 



w{h-^ji:^), (21) 



where A, A' , B' , jjgg and fitting parameters. Notice that in the above expression 

the inverse size (in impact parameter) figg is defined in the same way as in the original 
model, i.e., ^gg = y/fJ'qqfigg- In deriving the expressions (pUj) and (PT|) we have used the 
fact that the main contribution to the asymptotic behavior of hadron-hadron total cross 
sections comes from gluon-gluon semihard collisions, since g{x) ^ q{x) at small-x values. 
Therefore it is enough to build instrumental quark-quark and quark-gluon parametrizations 
for the expected high-energy behavior of the pp and pp amplitudes and to compute only 
the gluon-gluon contribution by means of the calculation procedure described in the last 
section. For example, the term ln(s/mp is naturally explained by the presence of a massive 
gluon in the qg — > qg subprocess. In this way the chosen eikonals reflect exactly the terms 
that come from such cross section calculations. We also have not considered the effect of 
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dynamically generated quark masses in the subprocesses involving quarks. This approach 
involves an extra parameter {rriq ^ 250 — 300 MeV), but we believe that its effect is smaller 
compared to the dynamical gluon mass one. 

The gluon-gluon contribution dominates at high energy and determines the asymptotic 
behavior of the total cross section. In our model we associate the gluon eikonal term Xgg{b, s) 
(see the expression (HH)) with the cross section a^^^{s): Xggipis) = a^^^ {s)W {b; figg) . 
Hence the gluon eikonal contribution includes gg —>■ gg subprocesses with colour nonsinglet 
exchange in all possible channels. The cross section o'^^^{s) is written as 

^fri') = C'[' drFgg{r)a^^^{s\ (22) 

where Fggij) is the convoluted structure function for pair gg, a^^'^{s) is the subprocess cross 
section given by expression and C is a fitting parameter. In the above expression we 
have introduced the energy threshold s > 4m^ for the final state gluons, assuming that these 
are screened gluons, in a procedure similar to the calculation of Ref. The structure 

function Fggij) is written as 

Fgg{r) = \g®g]{r) = ^ g{x) g (£j , (23) 

where g{x) is the gluon distribution function, usually adopted as 

9{x) = Ng L-^, (24) 

where J = 1 + e and Ng = — e)(5 — e)...(l — e). In this definition the^ term ~ 

l/x^^"^ simulates the effect of scaling violations in the small x behavior of g{x) [10]. In the 
Regge language the quantity J, that controls the asymptotic behavior of crtot{s), is the so 
called intercept of the Pomeron. In fact, neglecting the variation with of the asymptotic 
expression ((Tj), it is possible to show that 

hm drFggir) a^^^s) ~ f /^V • (25) 

Hence the total cross section behaves asymptotically as a Pomeron power law s"'^^, and a 
consistent value of J can be determined by fitting forward quantities data through a Regge 
pole model. Recently, by means of an extended Regge model, some authors have determined 
the bounds for the soft Pomeron intercept imposed by the accelerator and cosmic ray data 



currently available |24l. |25| . These results are consistent with a Pomeron intercept J = 1.085 
(specifically, J — 1 = 0.085 ± 0.006 in the case of constrained bounds 2^), and corroborate 
the choice of the Pomeron intercept value adopted is this work. 

We ensure the correct analyticity properties of our model amplitudes by substituting 
s se"""^/^ throughout Eqs. ((201), P^j) and (j221)- For simplicity, we will refer to our 
QCD-inspired model with a dynamical gluon mass simply as the DGM model. 

IV. RESULTS 

In all the fits performed in this paper we use a fitting procedure, where the value 
of Xmin is distributed as a distribution with N degrees of freedom (DOF). The fits to 
the experimental data sets are performed adopting an interval ~ Xmin corresponding, 
in the case of normal errors, to the projection of the hypersurface containing 90% of 
probability. In the case of the DGM model (8 fitting parameters) this corresponds to the 
interval ~ Xmin — 13.36 26j. To determine the optimum value for the dynamical gluon 
mass and extract the best phenomenological values of the DGM model parameters, we 
follow a two step process. First, we select specific input values for the dynamical gluon 
mass and carry out global fits to all high-energy forward pp and pp scattering data above 
^/s = 10 GeV and to the elastic differential scattering cross section for pp at ^/s = 1.8 
TeV. These forward data sets include the total cross section {<7tot), the ratio of the real to 
imaginary part of the forward scattering amplitude (p), and the nuclear slope in the forward 
direction {B). We use the data sets compiled and analyzed by the Particle Data Group P|, 
to which we add the new E811 data on af^^ and fP^ at ^/s = 1.8 TeV |27|. The statistic 
and systematic errors of the forward quantities have been added in quadrature. The input 
values of the rrig have been chosen to lie in the interval [300, 800] MeV, as suggested by the 
value mg = 500 ± 200 usually obtained in other calculations of strongly interacting processes 
(see section II). Although no physical argument ensures that the optimum value of rrig lies 
in the chosen input mass interval, our global fit results indicate a minimum value just about 
rrig ~ 400 MeV. These results are shown in Fig. 1, where a general dashed curve is added 
to guide the eye. Roughly, taking a 5% variation on the minimal x^/DOF value indicated 
by the general curve, it is possible to estimate a dynamical gluon mass rrig ^ mtfoo MeV. 
This result is totally compatible with the ones of Ref. rrig = 370 MeV. 
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Next, in order to determine the parameters of the DGM model, we set the value of the 
dynamical gluon mass to nig = 400 MeV (optimal value) and carry out a global fit only to 
all high-energy forward pp and pp scattering data above y/s = 10 GeV, not including the 
elastic differential scattering cross section da^^/dt at ^/s = 1.8 TeV. The values of the fitted 
parameters are given in Table 1. The x^/DOF for this global fit was 1.075 for 188 degrees of 
freedom. The results of the fits to atot, P and B for both pp and pp channels are displayed in 
Figs. 2, 3 and 4, respectively, together with the experimental data. Within this procedure, 
the Tevatron differential cross section, as well as the Tevatron-run II and the CERN LHC 
ones, can be predicted by the DGM model. These predictions are shown in Fig. 5. Table 
II contains predictions for the forward quantities at these energies, where the quoted errors 
are the statistical errors due to the errors in the fitted parameters. 



V. CONCLUSIONS 



In this paper we have investigated the influence of an infrared dynamical gluon mass 
scale in the calculation of pp and pp forward scattering quantities through a QCD-inspired 
eikonal model. By means of the dynamical perturbation theory (DPT), we have computed 
the tree level gg gg cross section taking into account the dynamical gluon mass, and have 
shown that the IR divergences associated with the gluon-gluon subprocess cross section are 
naturally regulated by this dynamical scale. In order to make a connection between the total 
subprocess cross section d'gg{s) and the forward pp and pp quantities, we have developed a 
QCD-inspired eikonal model where the onset of the dominance of gluons in the interaction 
of high energy hadrons is managed by the dynamical gluon mass scale. Using this formalism 
it was possible not only to reduce the number of parameters of the model, but also to give a 
consistent physical explanation for each one. For example, in some recent papers on QCD- 
inspired models 0, O], the two arbitrary constants mo and as were assumed to be equal 
0.6 GeV and 0.5, respectively; in our approach the IR value of running coupling constant 
is driven by the dynamical gluon mass, i.e., its IR behavior depends on the value of nig. 
This connection permit us to decrease the number of parameters required to describe the 
hadronic experimental data. 

By means of a global fit to the forward pp and pp scattering data and to da^^/ dt data 
at a/s = 1.8 TeV, we have determined the best phenomenological value of the dynamical 
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gluon mass, namely nig ~ 400_]^qo MeV. Interestingly enough, this value is of the same order 
of magnitude as the value nig ^ 500 ± 200 MeV, obtained in other calculations of strongly 
interacting processes. This result corroborates theoretical analysis taking into account the 
possibility of dynamical mass generation and show that, in principle, a dynamical nonper- 
turbative gluon propagator may be used in calculations as if it were a usual (derived from 
Feynman rules) gluon propagator. 

With the dynamical gluon mass set at rUg = 400 MeV, we have performed a global fit 
only to the forward pp and pp scattering data, in the same way as is usually performed in 
the former QCD-inspired models. Our model allows us to describe successfully the forward 
scattering quantities 0"^*, p and B, as well as to predict the pp differential cross section at 
^/s = 1.8 TeV in excellent agreement with the available experimental data. These results 
show that the DGM model is well suited for detailed predictions of the forward quantities 
to be measured at higher energies. In particular, for the total cross sections to be measured 
at Tevatron-run II and CERN-LHC energies, the model predicts the values atot = 75.7 ±5.4 
mb and atot = 102.9 ±7.1 mb, respectively. Our central LHC value prediction is close to the 
central one in the Ref. namely atot = 108 mb. This relatively small difference reflects 
the fact that the dynamical gluon mass and its associated coupling constant successfully 
replace the values of the ad hoc parameters of the former QCD inspired models. However, 
if the pp total cross section is measured at the LHC with a precision up to 5%, a selection 
between these QCD models may be possible. 

In summary, we argue that the QCD-inspired eikonal model with a dynamical IR mass 
scale provides an useful phenomenological tool to the study of the hadron-hadron diffractive 
scattering, where a purely perturbative QCD method is inadequate. 
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TABLE I: Values of the parameters of the DGM model resulting from the global fit to the forward 
pp and pp data. The dynamical gluon mass scale was set to rUg = 400 MeV. 



C (12.097±0.962)xl0-3 

Mgg [GeV] 0.7242±0.0172 

A 6.72±0.92 

Hgq [GeV] 1.0745±0.0405 

A' (4.491±0.179)xl0-3 

B' 1.08±0.14 

C- 3.17±0.35 

IJ,- [GeV] 0.6092±0.0884 



TABLE IL Predictions of the pp and pp forward scattering quantities crtot, p and B for the Fermilab 
Tevatron run-II (TEVII) and the CERN LHC energies. 

TEVII [L96 TeV] LHC [14 TeV] 
aZ, 4ot [mb] 75.7±5.4 102.9±7.1 
(jPP, fpP 0.129±0.009 0.114±0.005 
BPP, BPP [GeV-2] 16.97±0.99 19.36±1.12 
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FIG. 1: The /DOF as a function of dynamical gluon mass nig. 
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FIG. 2: Total cross section for pp (solid curve) and pp (dashed curve) scattering. 
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FIG. 3: Ratio of the real to imaginary part of the forward scattering amphtude for pp (sohd curve) 
and pp (dashed curve) scattering. 
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FIG. 5: Predictions for the elastic differential scattering cross sections at ^/s = 1.8,1.96 and 14 
TeV. In our model the channels pp and pp are not distinguished at high energies. The data points 
are from E710 |28|. 
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